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This paper is a review of heavy quarks in lattice gauge theory, focusing on methodology. It includes a status 
report on some of the calculations that are relevant to heavy-quark spectroscopy and to flavor physics. 



1. MOTIVATION 

The study of flavor- and CP- violation is a vi- 
tal part of particle physics [1]. Often lattice QCD 
is needed to connect experimental measurements 
to the fundamental couplings of quarks, which, in 
the Standard Model, are elements of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix. Usually the 
test of CKM is drawn as a set of constraints on 
the apex of the so-called unitarity triangle (UT). 
The Particle Data Group's version [2] is shown 
in Fig. 1. Apart from the wedge sin 2/3, theoreti- 
cal uncertainties dominate, and everyone wonders 
whether they are reducible and, if so, reliable. In- 
deed, as Martin Beneke [3] put it at Lattice 2001, 
the "Standard UT fit is now entirely in the hands 
of Lattice QCD (up to, perhaps, |V„b|)." 

The needed hadronic matrix elements are 
among the simplest in lattice QCD, so wc can 
hope to carry out a full and reliable error anal- 
ysis. Two criteria are key: First, there must be 
one stable (or very narrow) hadron in the initial 




Figure 1. Unitarity triangle from Ref. [2]. 



state and one or none in the final state; second, 
the chiral extrapolation must be under control. 
Such quantities can be called gold-plated, to re- 
mind us that they are the most robust. (They 
are conceptually and technically much simpler 
than non-lcptonic decays, or resonance masses 
and widths.) Moreover, realistic, unquenched 
simulations for gold-plated quantities now seem 
to be feasible [4,5]. 

Much is at stake. Fits to the CKM matrix are 
often described as over-constraining the Standard 
Model, but they are really a test. With (reliable) 
error bars of a few percent one can imagine a pic- 
ture like Fig. 2. The base, the angle 7, and the 
left side come from B (and other) decays that 
proceed at the tree level of the electroweak inter- 
actions. They can be considered to measure the 
UT's apex. The side labeled "V^t^ is obtained 
from the frequency of B°-B° mixing; the angle (3 
through the interference of decay with and with- 
out mixing. Standard mixing proceeds through 
box diagrams, so non-Standard processes could 
compete. Thus, this side tests the CKM theory. 
It could fail the test by missing the apex, thereby 
inspiring a grand revision of the program of flavor 
physics: "The CKM Matrix Reloaded." 
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Figure 2. A possible unitarity triangle in a few 
years. The scatter covers approximately 5<7. 
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Further motivation for lattice QCD with heavy 
quarks has emerged this year. The BaBar [6], 
CLEO [7] and Belle [8] experiments have ob- 
served the lowest-lying J p = + and 1 + states 
in the D s system. They lie slightly below the 
D(*>K thresholds, so they are narrow, decaying 
via isospin violation to D s ir. Several groups 
have started to calculate the masses with lattice 
QCD; see Sec. 6.2. The SELEX experiment re- 
ports not one or two but five baryons that pre- 
sumably contain two charmed quarks [9] . Corrob- 
oration from lattice QCD would be useful. A pilot 
study appeared during the conference [10]. 

The rest of this paper is organized as follows. 
In Sec. 2 the problem of heavy quarks in lattice 
gauge theory is deconstructed into three elements: 
a discretization of the heavy quark(s), an effec- 
tive field theory for describing cutoff effects, and 
a set of calculations of the short-distance behav- 
ior for matching to QCD. The most widely used 
methods are then cast into this framework. Sec. 3 
reviews the effective field theories needed to tame 
heavy-quark discretization effects, leading up to 
a semi-quantitative comparison of the methods. 
The issue of renormalon shadows is re-examined 
in Sec. 4. A new technique is discussed in Sec. 5. 
Sec. 6 discusses the mass spectrum and matrix el- 
ements needed for the Standard CKM fit. A short 
summary and outlook are in Sec. 7. 

2. METHODS FOR HEAVY QUARKS 

To connect lattice calculations of heavy-quark 
physics to experiment, we must confront three 
principal concerns: the quenched approximation, 
heavy-quark discretization effects, and the chiral 
extrapolation. The quenched approximation is 
the most worrisome to experimenters and phe- 
nomenologists, because its error is unquantifiable. 
Fortunately, it seems to be going away (at last). 
With 3 improved staggered sea quarks we find 
agreement with experiment for a wide variety 
of gold-plated masses and decay constants [4,5]. 
Some issues should be clarified [11], but it is 
clearly time to face the other two problems. 

Heavy-quark discretization effects arc vexing 
because, with available computing resources, the 
bottom quark's mass in lattice units is not small, 



mba ^ 1. (The charmed quark's mass is also not 
very small.) When ma <C 1, as for light quarks, 
one can control and quantify cutoff effects with 
a Symanzik effective field theory [12]. This kind 
of idea can be extended to heavy quarks [13,14]. 
The ensuing insights then allow us to compare 
and contrast the size and nature of heavy-quark 
discretization effects in the various techniques. 

The chiral extrapolation arises, as elsewhere in 
lattice QCD, because quark masses as small as 
those of the up and down quarks are not feasible 
computationally. The key issue [15] is whether 
the numerical data exhibit the curvature char- 
acteristic of chiral perturbation theory (%PT). If 
so, we can trust to guide the extrapolation. 
If not, the extrapolation is still necessary, but 
its uncertainty becomes difficult to estimate. In- 
deed, in at least one important example this un- 
certainty has been underestimated [16]. 

There are many ways to treat heavy quarks 
in lattice gauge theory. It is instructive to de- 
construct them into three essential elements, il- 
lustrated in Table 1. The first is the choice of 
discretization — a lattice fermion field and action. 
By "Wilson quarks" I mean four-component lat- 
tice fermion fields with any action using Wilson's 
solution of the doubling problem [17]. "Static 
quarks" are essentially one-component [18]; spin- 
dependent effects arc incorporated as inser- 
tions [19]. "Non-relativistic quarks" have two 
components [20], and the action is a discretiza- 
tion of a non-relativistic field theory [21]. Lat- 
tice gauge theories with any of these quark fields 
are compatible with the heavy-quark symmetries 
that emerge in QCD [22,23]. 

To run numerical simulations one needs only 
the discretization. To interpret the simulation 
data, however, one needs more. A typical simula- 
tion has unphysical values of some of the scales, 
such as lattice spacing a ^ and light quark mass 
much larger than the down quark's mass. A the- 
oretical framework is needed to get from simula- 
tions of lattice gauge theory (LGT) to the target 
quantum chromodynamics (QCD). A powerful set 
of frameworks is based on effective field theories. 
Indeed, most theoretical uncertainties in lattice 
QCD can (and should) be controlled and quanti- 
fied with effective field theories [24]. 
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Tabic 1 



The main elements for incorporating quarks with mass toq 3> Aqcd into lattice gauge theory. 



heavy-quark discretization 


effective field theory tools 


rcnormalization & matching 


(improved) Wilson quarks 

static quarks (+ insertions) 
non-relativistic quarks 


Symanzik LE£: for m,Qa <C 1 
for rriQa f£ 1 

HQET (for qQ) 
NRQCD (for QQ) 


perturbative: tadpole tree-level 

1- or 2-loop 
non-perturbative 
combination: Za = Zy P 



In the context of discretization effects, the best- 
known effective field theory is Symanzik's the- 
ory of cutoff effects [12]. The central concept 
is a local effective Lagrangian (LE£) of a con- 
tinuum field theory. It is also possible, when a 
quark is heavy, to use the (continuum) heavy- 
quark theories HQET (heavy-quark effective the- 
ory) and NRQCD (non-relativistic QCD) [13,14]. 
The next section discusses these ideas in detail. 
Here I note only that I distinguish HQET from 
NRQCD by the way they classify interactions, as 
dictated by the physics of heavy-light and heavy- 
heavy hadrons, respectively 

Effective field theories separate dynamics at 
short distances (a, 1/mq) from those at long dis- 
tances (1/Aqcd, 1/?ti<3, L). To rcnormalizc LGT 
(or to "match" LGT to QCD), one must calcu- 
late the short-distance mismatch between LGT 
and QCD, and tune it away. It is natural to cal- 
culate the short-distance behavior with perturba- 
tive QCD, owing to asymptotic freedom, and au- 
tomated multi-loop calculations for improved ac- 
tions arc feasible [25]. Non-perturbative methods 
for matching are also available [26,27], although 
most of the results apply only when rriQa <c 1. 
Combinations, in which some (or most) of the 
rcnormalization is computed non-pcrturbatively 
[28-30] are possible too. For example, the nor- 
malization factor Zy(m,Qa) for the vector cur- 
rent is easy to compute non-perturbatively for all 
raqa. Then pa = Za/Zv is tadpole-improved, 
and explicit examples [14] exhibit very small one- 
loop corrections and mild mass dependence. 

With these ideas in mind, we can summarize 
quickly the most common methods for lattice B 
physics, in each case working across Table 1. 

The first method [31,32] dates back to the 1987 
Scillac conference. It starts with Wilson quarks, 
nowadays almost always with the improved clover 
action [33]. Cutoff effects are described through 



the original Symanzik LE£ with rriQa <C 1. The 
advantage is that non-pcrturbativc matching cal- 
culations derived for light quarks [27] may be 
taken over. The disadvantage is that the heavy- 
quark mass in the computer is artificially small 
(to keep mga < 1), so matrix elements of heavy- 
light hadrons are extrapolated in 1/mq from 
rtiQ ~ m c up to 77i &. This method does not work 
for heavy-heavy systems, because toq dependence 
in these systems, from NRQCD, is not character- 
ized by powers of I/thq. 

There are several names for this method in 
the literature. Adherents prefer "the relativistic 
method" or, simply, "QCD." The first name over- 
looks the major discretization effect, which stems 
from violations of relativistic invariancc. The sec- 
ond is not helpful: all methods start with lattice 
gauge theory, and use an effective field theory to 
define the matching and to interpret cutoff effects. 
Especially after the extrapolation in I/tuq, this 
technique is not closer than the others to contin- 
uum QCD. In the rest of this paper, I call this 
method "the extrapolation method." 

At Seillac, Eichten introduced the method now 
known as lattice HQET [18]. Indeed, this work 
spawned the continuum HQET. It starts with a 
discretization of the static approximation. Then 
1/rriQ corrections are treated as insertions [19]. 
The Symanzik LE£ is then nothing but the con- 
tinuum HQET [34]. Non-perturbative renormal- 
ization methods have been worked out in the 
static limit [35,36] but are sketchy for 1/toq in- 
sertions. The static approximation and, hence, 
lattice HQET do not apply to quarkonium. 

Another technique, also introduced at Scillac, 
is called lattice NRQCD [20]. It first derives 
a non-relativistic effective Lagrangian for heavy 
quarks in the continuum, and then discretizes. 
The deviation from QCD, both from the non- 
relativistic expansion and from the discretization, 
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can be described with continuum heavy-quark 
theory: HQET for heavy-light or NRQCD for 
heavy-heavy systems. For the latter several inter- 
actions are needed (for few-percent accuracy), so 
the lattice action can be rather complicated [37]. 
For this reason most of the matching calculations 
are perturbativc. The same lattice action is used 
for heavy-light and heavy-heavy systems, so the 
latter can be used to test its accuracy, before pro- 
ceeding to B and D physics. 

The Fermilab method [38] is a synthesis. Like 
the extrapolation method, it starts with (im- 
proved) Wilson fermions. It confronts heavy- 
quark discretization effects in two ways. One is 
to extend the Symanzik effective field theory to 
the regime toqo 1. The other is to use (contin- 
uum) HQET and NRQCD to describe both short 
distances, a and 1/tbq. This is possible because 
Wilson fermions have the right heavy-quark sym- 
metries. Renormalization and matching is usually 
partly non-perturbative and partly perturbative. 
Like lattice NRQCD, the Fermilab method can be 
used for heavy-light and for heavy- heavy systems. 

Heavy-light hadrons also contain light valence 
quarks. In the past, Wilson quarks were usu- 
ally used, but Wingate el al. have shown the ad- 
vantages of naive quarks [39] . Several difficulties 
associated with doublers are absent because the 
naive quark is bound to a heavy quark. Naive 
quark propagators are obtained from staggered 
quark propagators by undoing the spin diagonal- 
ization. With staggered quarks one can reach 
much smaller light quark masses (m q ~ 0.2m s ) 
than with Wilson quarks (m q ~ 0.5m s ) [5,11]. 

The main issue with the light quarks is the chi- 
ral extrapolation. Heavy-meson chiral perturba- 
tion theory describes the dependence of heavy- 
light masses and matrix elements on the light- 
quark mass [40]. It can be extended to par- 
tially quenched simulations [41]. For quarkonium 
the chiral extrapolation is less well developed. 
The virtual process most sensitive to the light- 
quark mass is dissociation, which could have a 
strong effect on states just below threshold, such 
as ip' ~ * DD — > ip'. This process is completely 
omitted in the quenched approximation. 

To close this section, let me list the (perceived) 
problems with the four main treatments of heavy 



quarks. The extrapolation method must balance 
the contradictory requirements of keeping both 
Aqcd/it-q and rriQa small. In practice it is not 
clear that both sources of uncertainty, and their 
interplay, are controlled. In lattice HQET or lat- 
tice NRQCD, power-law divergences arise from 
the short-distance behavior of higher-dimension 
interactions and remain at some level when per- 
turbative matching is used. There are also the 
usual uncertainties from truncating perturbation 
theory. The Fermilab method, being based on 
Wilson fermions, does not have power-law diver- 
gences, but the perturbative part of the matching 
is said to suffer from a remnant called "renor- 
malon shadows." Another problem is that the 
a dependence, though smooth, is not easily de- 
scribed in the regime rriQa ~ 1. These problems 
are assessed at the end of the next section, after 
discussing the theory of cutoff effects in detail. 

3. CUTOFF EFFECTS 

In this section I would like to discuss heavy- 
quark discretization effects in more detail. To do 
so, I shall use effective field theories to separate 
short- and long-distance dynamics. Discretiza- 
tion effects are at short distances, and here the 
methods differ. One can then assess the uncer- 
tainties from the short-distance behavior, without 
precisely knowing the long-distance effects. 

Although effective field theories appear in 
many contexts, it is perhaps worth recalling some 
of the basics. At energies A below some scale /z, 
particles with E > p, have small effects, sup- 
pressed by (A/E) n . The Coleman-Norton the- 
orem [42] shows that analytic properties of Green 
functions arc impervious to off-shell particles. 
This is because singularities appear only where 
particles go on shell. Suppose a heavy quark zig- 
zags into an anti-quark (emitting a hard quan- 
tum), and then emits a soft quantum before turn- 
ing back into a quark (absorbing a hard quan- 
tum). The anti-quark is far off shell. The ana- 
lytic structure is retained even if the anti-quark 
and hard propagators are reduced to a point. 
The resulting reduced diagram already looks like 
an effective-field-theory diagram. New fields can 
be introduced to describe the remaining particles 
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and. hence, their non-analyticity. To compensate 
for the omitted, off-shell particles, vertices require 
generalized couplings. This framework suffices, 
because field theory gives a complete description 
respecting analyticity, unitarity, and the underly- 
ing symmetries [43] . 

Asymptotic freedom provides another way to 
establish effective field theories for underlying 
gauge theories, namely to all orders in the gauge 
coupling. In this way we have a rigorous proof 
(for Wilson quarks) of Symanzik's theory of cutoff 
effects [44], and a (less rigorous) proof of the static 
field theory on which HQET is based [45] . Given 
the general arguments [42,43] sketched above, 
it is hard to see where these would fail non- 
pcrturbatively. In particular, it seems implausi- 
ble that confinement would be different with two- 
component heavy-quark fields instead of four- 
component Dirac fields. 

With heavy quarks vtlq 3> Aqcd, so zig-zags 
and pair production are suppressed in bound 
states. Two-component fields hi + ^ (hi ') suffice 
to describe quarks (anti-quarks). Here v is a la- 
bel, denoting a four-velocity close to that of the 
hadron containing the heavy quark(s). In systems 
with one heavy quark, the heavy quark's velocity 
deviates from v by an amount of order Aqcd /tuq. 
The field theory describing this physics with hv 
fields is HQET. If there are two heavy quarks, 
one has a binary system. The two heavy quarks 
rotate slowly about their center of mass, generat- 
ing long-distance scales of order mqv", where v ~ 
a s (rn>o) denotes the relative three- velocity. In the 
bb (cc) system v 2 w 0.1 (0.3) [37]. Now the field 
theory with hi ±] fields is NRQCD. 

HQET and NRQCD are based on an equiva- 
lence between QCD and an effective Lagrangian. 
One may write 



C 



QCD 



Cm 



ght 



-HQ, 



(1) 



where the symbol = can be read "has the 
same matrix elements as." We cannot use = 
("equals"), because zCqcd has Dirac fields for 
heavy quarks and £hq has heavy-quark fields 
h^\ The effective Lagrangian 



Cm = ^2Ci{m Q ,m Q / 'n)Gi(n), 



(2) 



where the Cj are couplings, also known as short- 
distance coefficients. The operators Oi describe 
long-distance dynamics and, thus, bring in soft 
scales Aqcd, WqF, etc. The renormalization 
scale of the effective theory is /x, chosen so that 
(Aqcd, itt-qv) < fi < rtiQ. The difference be- 
tween HQET and NRQCD reflects the physical 
differences between heavy-light and heavy-heavy 
systems, classifying the terms in Eq. (2) accord- 
ing to powers of the respective small parameters, 
Aqcd/toq and (m Q v)/mQ = v. 

It is easiest to see the difference between the 
two by looking at the leading and next-to-leading 
terms. Let us choose the frame v = (1,0) and, for 
brevity, suppress the velocity label. For HQET 
one counts powers of Aqcd / m Q , 



4jq = -M m i 



C 



(i) 

HQ 

(2) 
HQ 



hD 2 h 



D 4 )h, 
hill ■ Bh 



2m 2 
hD ■ Eh 

ZD- 



2m,2 



8m\ 



+ z s 



hiT, ■ [D x E]h 



(3) 
(4) 

(5) 



and in general £hq consists of all interactions of 
dimension s+4. The Zi are dimcnsionlcss versions 
of the Ci\ they depend on uiq/ii with anomalous 
dimensions. Thus, HQET has both power-law 
and logarithmic ttiq dependence. 

For NRQCD, on the other hand, one counts 
powers of v. As explained in Ref. [37], D ~ itiqv, 



B ~ m 2 v 4 , E 



4jq = -M m i 



Av 3 . Then 



D 4 )h 



,(2) 
'HQ 



zb- 



hiH ■ Bh 



hD 2 h 
hD ■ Eh 



zu- 



2m 2 
hiH ■ [D x E]h 



8m| 



+ 

h (D 2 ) 2 h 
8ml 



(6) 
(7) 



In general consists of all terms that scale 

like v s+2 . Eqs. (6) and (7) contain the same terms 
as Eqs. (3)-(5), etc., but rearranged. For exam- 
ple, the kinetic energy is an essential effect in 
NRQCD, but a non-leading correction in HQET. 

The rest mass m\ and kinetic mass ttt-2 should 
be thought of as short-distance coefficients. 
When describing QCD with £hq, mi = mi- If 



6 



the operators are renormalized in the MS scheme, 
this mass is (in perturbation theory) the pole 
mass. When we turn to LGT below, mi =/= m%. 

I believe that the key to reliable calculations 
lies in firmly incorporating these heavy-quark 
ideas into LGT. After all, almost nothing is 
known about heavy quarks (in bound states) 
without these and other scale-separation tech- 
niques. In particular, the machinery outlined here 
can be applied to any underlying field theory with 
heavy-quark symmetry, including LGT with Wil- 
son, static, or non-relativistic quarks. Many read- 
ers will agree with such a strategy, but others 
seem to believe that effective theories should be 
avoided. If so, they overlook the fact that the al- 
ternative theory of cutoff effects — Symanzik's — is 
(!) an effective field theory. 

Let us turn, then, to the Symanzik theory, with 
an eye to nuances that arise for heavy (Wilson) 
quarks. To describe cutoff effects, Symanzik [12] 
introduces a local effective Lagrangian (LE£). 
We can write £lgt = £sym, and 

£ Sym = ^ trL^FH - £ q f (p + my)?/ (8) 

/ 



and soft terms [38,46,24]. One finds 



2.9 2 



i 

The first line gives continuum QCD, with renor- 
malized coupling g 2 and quark masses my. The 
sum in Eq. (8) accounts for cutoff effects. Quarks, 
even heavy ones, arc described by fields qy satis- 
fying the Dirac equation. The renormalization 
scale fi separates the short distance a from long 
distances, principally 1/Aqcd- For a light quark, 
1 /m y is another long distance, so it is sensible and 
accurate to expand the coefficients i-Q in (mya)™. 
In this way one obtains an expansion in a that is, 
however, a consequence of separating scales when 
the only short distance is a. 

For a heavy quark Q, on the other hand, 1/toq 
is (compared to 1/Aqcd) a short distance. Hence, 
it does not make sense to expand the Ki in toqo. 
The expansion is also inaccurate when rriQa 1. 
Instead of expanding, one should rearrange the 
terms in Eq. (8) to collect the largest terms from 
the sum. They are of the form QXlp^D^Q + 
h.c, which can be eliminated with the equations 
of motion in favor of explicit mass dependence 



£sym = C-F.q - Q I mi + 7 4 £>4 

+ small corrections, 
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(9) 



The incorrect normalization of the spatial kinetic 
energy captures the breaking of rclativistic invari- 
ance when m,Qa ^ 1. This modified LE£ is on 
the same footing as Eq. (8). Thus, one sees that 
for toqo 1 it is neither lattice gauge theory nor 
the Symanzik LE£ that breaks down. Instead, 
the split "QCD + small corrections" suggested in 
Eq. (8) is lost in Eq. (9). 

The apparent obstacle can be overcome by us- 
ing HQET and NRQCD to separate the short dis- 
tances (a, 1/toq) from (Aqcd, m.QV n ). The logic 
and structure for LGT is the same as for QCD 
[13]. With two short distances, the Wilson coef- 
ficients depend on the dimensionless ratio rriQa: 

Ci(m Q ;n/m Q ) -> Ci(mQ,m Q a;{cj}; n/mq). (10) 

As shown, they also depend on improvement cou- 
plings Cj in the lattice action, as do the coeffi- 
cients Ki in Eq. (8). 

The rest mass mi and kinetic mass mi in 
Eq. (9) are the same as Eqs. (3)-(6). But in 
HQET and NRQCD m\ multiplies the conserved 
number operator hh: it contributes additively to 
the mass spectrum but not at all to matrix ele- 
ments. So, although it is possible to introduce a 
new parameter to tunc iiij = m 2 [38] , it is also un- 
necessary (assuming tuq 3> Aq CD ). In this way, 
the heavy-quark theories make sense out of the 
short-distance behavior of Wilson quarks, where 
the Symanzik theory does not. 

With HQET and NRQCD descriptions of 
static, non-relativistic, and Wilson quarks, we 
are in a position to make comparisons. Heavy- 
quark expansions for QCD and LGT can be devel- 
oped side- by-side [47], showing explicitly [14] that 
heavy-quark discretization effects are lumped into 

(11) 



QCD 



^•=c^(te})-c; 

Solving the equations 5Ci = for the couplings Cj 
yields on-shell improvement conditions. For cur- 
rents there is an analogous description leading to 
similar mismatches and to normalization factors 

Zj = C? CD /CJ GT , 



(12) 
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where Cy is the matching factor from the under- 
lying theory to the effective theory. The resulting 
expressions for the cj and Zj no longer depend on 
the renormalization scheme of £hq: they match 
LGT directly to QCD. 

To turn these results into semi-quantitative es- 
timates, one can look at the mismatches in the 
coefficients. Let us focus on the action for illus- 
tration. In the extrapolation method the quark 
mass is identified with the rest mass, niQ = m\, 
so the kinetic energy introduces an error 



SC. 



(extrapolation) 



(O2) 



p 



p 



2ui2 2m\ 



(13) 



iQ=mi 



where p is a soft scale. The chromomagnetic mis- 
match SCb(Ob) is similar. In lattice NRQCD the 
quark mass is adjusted non-perturbatively to the 
kinetic energy: mg = m^. The first error is in 
spin-dependent effects. With /-loop matching 



j£ (NRQCD) 



(Ob) 



I 



1 



P 



2mq 



(14) 



has the right asymptotics in niQa. The power-law 
divergence comes from the short-distance part 
of higher-dimension spin-dependent interactions. 
The Fcrmilab method chooses tuq = 7712 and, so, 
eliminates the error (13). The leading error is 



SC 



(Fermilab) 
B 



(Ob) 



a l s +1 pa 



2(1 + moa) 



(15) 



rnQ=m 2 



again with /-loop matching of the chromomag- 
netic energy. The existence of a continuum limit 
controls the a — > limit. 

The error estimates of Eqs. (13)— (15) are plot- 
ted in Fig. 3, taking p = 700 McV (~ A QC d), 
m c = 1400 McV, m b = 4200 McV, a s = 0.25, and 
1 = 1. The vertical lines indicate the lattice spac- 
ings of the unquenched MILC ensembles [5,48]. 
The results substantiate some of the perceived 
problems. For both charmed quarks (gray curves) 
and for bottom quarks (black curves), the ex- 
trapolation method has largest discretization ef- 
fects, 10-20% at feasible lattice spacing. Lattice 
NRQCD (with perturbative matching) has an in- 
creasing error as a — > 0, as noted in the first paper 
on NRQCD [20]. The Fermilab method's heavy- 
quark discretization effects are not a pure power 
of a at feasible lattice spacings, where toqo ~ 1. 



Calculations of matrix elements have similar er- 
rors from matching and further errors from the 
normalization factors. With lattice NRQCD and 
the Fermilab method the latter are, till now, of or- 
der a 2 s from one-loop matching. Two-loop match- 
ing is underway [25]. Matrix elements from the 
extrapolation method typically have leading nor- 
malization errors of order a s m,Qa . Some of the 
"improvement" contributions, taken over from 
light quarks, diverge in the heavy-quark limit. 
They are worrisome and should be avoided. 

Fig. 3 reveals two other interesting features. 
First, the extrapolation method has smaller dis- 
cretization effects than lattice NRQCD only at 
unfeasibly small lattice spacings. Second, dis- 
cretization errors from the light quarks, here 
taken to be i(pa) 2 , seem to dominate in lattice 
NRQCD and the Fermilab method. That sug- 
gests that one could carry out a continuum ex- 
trapolation to get rid of these effects, tolerating 
a percent-level bias from the heavy quarks. 

One can take care of the discretization effects 
of the extrapolation method by taking the contin- 
uum limit first [49], but it is dangerous. On the 
lattices available today (especially in unquenched 
simulations), one is driven to artificially small val- 
ues of the heavy-quark mass. At some point, 
however, the 1/toq expansion breaks down. The 
threshold effects that allow pair production and 
zig-zags are smeared out by the momentum dis- 
tribution of gluons inside the hadron, so, unfor- 
tunately, the breakdown will not be obvious. 




a (fm) 



Figure 3. Comparison of errors. Black: Q = b: 
gray: Q = c. Solid: extrapolation method; 
dashed-dotted: lattice NRQCD; dashed: Fermi- 
lab method; dotted: gluons and light quarks. 
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The above discussion skips over lattice HQET. 
With perturbativc matching, the kinetic, as well 
as the chromomagnctic, error is like (14). For 
this and signal-to-noise reasons, this method has 
lately fallen from favor. The latter problem now 
seems to be solved [50], and there is progress on 
non-pcrturbativc matching [35]. The way I have 
set up HQET makes it clear, however, that non- 
perturbative matching, say in a volume of finite 
size L, should work with minor modifications for 
static, non-relativistic, and Wilson quarks. Non- 
perturbative matching in finite volume will intro- 
duce uncertainties of order Aqcb/itt-qL. To con- 
trol these L cannot be too small. 

4. RENORMALON SHADOWS 

We now turn to renormalon shadows [51], 
which are perceived to be a problem in the Fer- 
milab method. 

Renormalons are power-law ambiguities that 
arise in large orders of perturbation theory, when 
mass-independent renormalization schemes are 
used. Several years ago, Martinelli and Sachrajda 
considered the problem of power corrections [52]. 
They concluded that a sluggish cancellation of 
renormalon effects could make it difficult to take 
power corrections into account. Frankly, I am 
suspicious of sweeping conclusions that hinge on 
renormalons, because in other renormalization 
schemes they do not arise [53]. Let us see, there- 
fore, whether it is possible to arrive at similar 
conclusions without reference to renormalons. 

Suppose one can measure experimentally (or 
compute non-perturbatively) V and Q. Both are 
described in an effective field theory including a 
power correction: 

V(Q) = CriQ/riiOx) + B v {Q/^)(0 2 )/Q, (16) 

Q(Q) = C Q (Q/n){Ox) + B Q (Q/n)(0 3 )/Q, (17) 

where the Cs and Bs are short-distance coeffi- 
cients. Q is a hard physical scale, and /z is the sep- 
aration scale. By assumption, the same effective- 
theory matrix elements (0\) and (O2) appear in 
V and Q. Now one would like to predict 

n(Q) = Cn(Q/n)(Ox) + B n (Q/(,)(0 2 )/Q, (18) 

given V, Q, and approximate short-distance co- 
efficients. What is the uncertainty in 72,? 



One could solve Eqs. (16) and (17) for the (O,), 
but renormalons could enter through the renor- 
malization scheme for the effective field theory. 
Instead, just use linear algebra to obtain: 



C 



■R 



v_ 



Q_ 



B n [T/C v - Q/C a ] 



(19) 



B-p/C-p - Bq/Cq 

where B^ is a combination of the Bs and Cs. The 
second term is the power correction: the combina- 
tion V/C-p — Q/Cq is formally of order Aqcd/Q- 

To estimate the uncertainty in 1Z, let us as- 
sume that the Cs have been calculated through / 
loops and the Bs through k loops. Then the cus- 
tomary assumption is that the uncertainty in the 
(leading-twist) first bracket is 0(a l s +1 ) and in the 
power correction 0(a k+1 Aqcd/Q). Moreover, it 
seems worthwhile to include the power correction 
as soon as Aqcd/Q > a£ +1 . 

But the leading-twist parts of V/Cp — Q/Cq 
do not cancel perfectly when Cp and Cq are cal- 
culated perturbatively. There is a shadow [51] of 
order a l s +1 , commensurate with the uncertainty in 
the leading-twist term. Unless a4 +1 <C Aqcd/Q, 
the shadow obscures the desired power correction, 
which is especially likely if V and Q are dissimi- 
lar [52]. Then the second term in Eq. (19) would 
not improve the accuracy of TZ. 

To see if this reasoning applies to any of the 
heavy-quark methods, one must check if the same 
kind of cancellation is needed. In lattice HQET, 
where 1/mg corrections are insertions, the fore- 
going analysis goes through without change, so 
the shadow is an issue. The Fermilab method ap- 
plies HQET in a different way, so the conclusions 
can be (and are) different. The HQET descrip- 
tion gives formulas of the form [13,14] 

$ LGT = C LGT [$Qo + B LGT^ /mQ ^ _ (2Q) 
^QCD = C QCD + B QCD^ /mQ ^ (21) 

We want to know the uncertainty, when Z = 
C^ /C LGT is computed through I loops, and 
SB = B LGT - £?Q CD is matched through k loops. 
One finds the relative error 

1 _ _ = " ^22., (22) 



$QCD 



z 



where Z^ is the Z-loop approximation to Z , and 
§Z( l+1 ) = Z - ZW. The first term has truncation 
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error a z s +1 , the second aJ +1 AQ CD /mg. There is 
no shadow contribution, because the power cor- 
rection is not obtained by explicit subtraction. 
Instead, it is present all along, and LGT is ad- 
justed to hit its target, QCD. 

One can always worry that the first uncalcu- 
lated coefficient in Z^ lJrV > or J£>(' +1 ) is large. But 
those are the errors discussed in Sec. 3. They 
have nothing to do with renormalons or shadows. 

5. A NEW TECHNIQUE 

Among this year's several methodological de- 
velopments, I would like to discuss a new finite- 
volume technique [54,55]. The idea is to calculate 
a heavy-light observable $ in a sequence of finite 
volumes of size Lq, 2Lq, .... Then 

$(<x>) = $(L )a(L )a(2L ) ■ ■ ■ , (23) 

where a(L) = $(2i)/$(L). To see how it works, 
let us follow the published work [54,55]. With 
Lq = 0.4 fm the box is small enough so that 
rriba <C 1 is feasible (though not cheap!). Thus, 
$(Lo) can be obtained in the continuum limit. 
One obtains a(2^~ 1 Lq) by taking the continuum 
limit first, and then extrapolating in 1/toq from 
2 3 /nib to 1/mb- "Infinite" volume is reached after 
two steps, where L = 4Lq = 1.6 fm. 

Systematic uncertainties stem from the 1/toq 
extrapolations. Let us focus on meson masses, 
which, in HQET, are written M(L) = to + A(L), 
where to is the quark mass. Renormalization- 
scheme dependence cancels between to and A. A 
encodes the long-distance physics, so it depends 
on the box size L, but to does not. Thus, 

cr(L) = 1 + [A(2L) - A(L)] /to. (24) 

Because of the extrapolation, one really has 

( L o) = 1 + 7^J7 — >( 25 ) 

where "2 : '" expresses the uncertainty arising from 
the extrapolation. Because a is computed in the 
continuum limit, these uncertainties arise from 
higher orders in 1/to (for j = 1) and very possibly 
from a breakdown of the heavy-quark expansion 
(for j = 2). When extrapolation is most danger- 
ous (larger j), the difference A(2 j L ) - A(2 j ^ 1 L ) 
nearly vanishes, according to the usual asymp- 



totic L dependence of hadron masses. 

The heavy-light decay constants can be ana- 
lyzed similarly, but HQET anomalous dimensions 
make the 1/to extrapolation less clean. The can- 
cellation mechanism still works, however. 

For bottomonium the extrapolations do not 
work as above. NRQCD says M (L) = 2m+B(L), 
where the binding energy B ~ tov 2 , with rela- 
tive velocity v ~ a s (to) . Extrapolations linear or 
quadratic in 1/to are, thus, not motivated. 

6. SPECTRUM AND CKM RESULTS 

This section surveys physics results, although 
less thoroughly than in the past two years [56,57]. 
It now makes sense to focus on unquenched cal- 
culations, which are just beginning. Here I focus 
on the quarkonium spectrum, the new D s states, 
semi-lcptonic decays, and B-B mixing. 

6.1. Quarkonium 

Fig. 4 shows the gross and fine structure of bot- 
tomonium [58], comparing quenched results to 
those on the MILC ensembles [5,48] with 2 + 1 
flavors. The lattice spacing is fixed by the spin- 
averaged IP- IS splitting, and the overall mass 
by the T(1S) level. In both parts of Fig. 4 one 
sees much better agreement with experiment for 
2 + 1 sea quarks. Note particularly the improve- 
ment in the excited levels 25, 35, and 2P, and 
the improvement in the Xbj splittings. 

We have similar results charmonium too. In 
this case, we can compare the hyperfine splitting 
toj/^— m Vc to experiment, finding an 18±2% de- 
viation [59]. This is better than in the quenched 
approximation, but still unsatisfactory. Note, 
however, that the heavy-quark chromomagnctic 
interaction is matched only at the tree level. One- 
loop matching [60,61] should fix the problem. 

The unquenched quarkonium spectrum allows 
us to fix the scales ro and r\ from first principles. 
Calculating the heavy-quark potential and fixing 
lattice units from the T spectrum, the MILC and 
HPQCD collaborations find [62] 

r = 0.46(1) fm, n = 0.32(1) fm, (26) 

at both lattice spacings, | and jj fm. Eq. (26) 
is the first QCD determination of ro; the value 
ro = 0.5 fm [63] comes from potential models. 
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6.2. D a (0+) and D*(l+) 

Now let us turn to the new positive-parity D s 
states [6-8]. The experiments find [7] 



m Ds (0 + ) -m Ds {0-) 

TO -Dj(l + ) _ W£)j(l~) 



350(1) MeV, 
351(2) MeV. 



(27) 
(28) 



These observations have inspired several lattice 
calculations [64-66]. The UKQCD Collaboration 
includes a simulation with n/ = 2, m™ a « m s , at 
fixed lattice spacing, a f» 1/10 fm [65]. Ref. [65] 
converts from lattice units to physical units with 
r = 0.55 fm. With r = 0.46 fm, from Eq. (26), 
their unquenched results become [67] 

m Ds (0+) -m Da = 494(28) MeV, (29) 
m Ds (l + ) -m Da {\-) = 445(28) MeV, (30) 

while m D; (r) - m B ,(0") = 140(3) MeV agrees 
with experiment. (Errors from non-zero a and 
imperfect unquenching are not reported here.) 

These splittings are too high, but in this case 
that is what one should expect. In Nature the 
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Figure 4. Bottomonium spectrum [58]. Open 
(closed) symbols show quenched (2 + 1) results, 
(a) Gross structure with lattice spacings a = g fm 
(middle) and a = jj fm (right). (b) Spin- 
dependent splittings with a 
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threshold D^if lies just above Z?i* ) (J+), so the 
virtual process £> S (J+) -> D^(J~)K -> D S (J+) 
pushes m £) (,)(J + ) down. UKQCD's sea-quark 
mass m™ a is slightly higher than the physical 
strange quark mass, so the threshold in simula- 
tion is higher, diluting its importance. If m s q ca 
would be reduced (eventually to raj), the thresh- 
old effect would reappear. Some evidence for 
this mechanism comes from the MILC ensem- 
bles [5,48] with lighter sea quarks, where we pre- 
liminarily do find lower splittings [66]. 

6.3. Semi-leptonic Decays 

Semi-leptonic decays determine the top two 
rows of the CKM matrix, because they are un- 
likely to be sensitive to non-Standard physics. 
Lattice QCD calculates the hadronic form factors, 
/+(g 2 ), f+(E n ), J-(w), etc., from matrix elements 
\k\V*\K), {-k\V»\B), (D*\J"\B), etc. 

To determine \V c b\ through the exclusive decay 
B — > D*lv : one needs the zero-recoil form factor 
Exact heavy-quark spin and flavor symme- 
tries would imply !F{\) = 1 [23]. The (quenched) 
state of the art is [30] 

lY-n — n niq+0.024 i n mfi +0. 003+0.000+0. 006 (o-t\ 
— U.»J.d_ 017 it U.UlO_ . 014 _o .016-0 .014) 

where the error bars are from statistics, HQET 
matching, non-zero lattice-spacing, chiral extrap- 
olation, and quenching. The total uncertainty 
of 4% hinges on the HQET description of LGT, 
so all uncertainties scale with .F(l) — 1, not ^(1). 

How will these uncertainties change in future 
unquenched calculations? The HQET match- 
ing error bar ±0.016 is reducible through higher- 
order improvement [68]. Of the rest, the chiral 
error bar ioloie ^ s the most subtle. The D* in 
the final state is unstable and, hence, not clearly 
gold-plated. In Nature, however, mr>* — hid is 
only slightly larger than m^, so the D* is nar- 
row, Td* ps 100 kcV [69]. In any feasible simu- 
lation, the artificially high mass of the light va- 
lence quark puts the D* above threshold. This 
effect has been worked out in xPT [70,41] and, in 
fact, leads to the io'o?6- With more experience 
and, hence, confidence in chiral extrapolations, 
one may hope to reduce (or at least symmetrize) 
this uncertainty. This would then be a 1% uncer- 
tainty on \Vcb\- 



11 



Charmless scmi-lcptonic decays can be used to 
determine \V ub \. The CLEO Collaboration [71] 
finds that B — * nli> is favored experimentally over 
B — > plv. That is fortunate for lattice QCD, 
because the 7r is stable and the p is not. Sev- 
eral unquenched calculations are in progress [72] . 
The final-state pion can reach high energy, E„ < 
2.6 GeV, which is another source of discretiza- 
tion errors. Ideas like moving NRQCD [73,74] 
promise to manage it. As usual one must worry 
about the chiral extrapolation. A nice LGT- 
oriented guide [75] for low comes from (par- 
tially quenched) xPT. One would like to treat 
energetic final-state pions as well, which, as far 
as I know, is an open problem. 

6.4. f B and B° q -B° q Mixing 

Neutral Bq-B® mixing plays a central role in 
B physics. (Here q labels the spectator d or s 
quark). The theoretical expression for the mix- 
ing frequency can be related to physics at the 
electroweak scale (and beyond), if one knows a 
QCD quantity called r)B-M q . Here 

M q = - 7 5 MMl - 7 5 )<?K), (32) 

and r\B is a short-distance coefficient. A4 q and tjb 
depend on the scheme for integrating out W , Z , 
and t, but rjsMq does not. The renormalization- 
group invariant value of the short-distance factor 
is t)b = 0.55. 

For historical reasons one usually writes 

B Bq = 3M q /8m 2 B J Bq , (33) 

f Bq = (0\h,l 5 q\B° q )/m Bq . (34) 

The separation into the decay constant fs and 
bag factor Bs q turns out to be useful. 

Am<2 is precisely measured [2], and Am s will 
be soon [76] . It is tempting to test the CKM ma- 
trix with Am s / Arrid- Then lattice QCD should 
provide the "SU(3)-breaking" ratios Ai s /Md or 

£ = £/£b, t f = f B Jf Bd , e B =B B JB Bi . (35) 

For several years, the conventional wisdom held 
that £ (and £/) would enjoy large cancellations in 
the uncertainties. Several uncertainties do cancel 
in the ratios, but one of the largest — from the 
chiral extrapolation of £/ — does not. 



A correct way to carry out the chiral extrapola- 
tion is, of course, to use xPT. In the case at hand, 
xPT provides the description f Bq = F[l + Af q ], 
Bb = B[l + AB q ], where F and B are constants 
and Af q and AB q are calculated from loop dia- 
grams. For the Bd meson one finds [77] 

Af d = -G+ [f IK) + \l{m 2 K ) + ±I(ml)] 

+ m 2 K f 1 + \ml[f 1 + fo], (36) 

ABd = —G- [\l(ml) + ±/(m 2 )] + m 2 K B x 

+ %ml[B 1 + B 2 ], (37) 

where G± = (l±3.g 2 )/(47r/) 2 , / = 130 McV, and 
g 2 is the B-B*-n coupling. x?T for the B s meson 
does not yield pion contributions /(m 2 ). A rough 
guide for g 2 comes from the D* width. Below I 
shall take g 2 = 0.35±0.15 to reflect the difference 
between the D and B systems. 

The dynamics of the QCD scale is lumped into 
the constants /; and Bi . The long-distance radia- 
tion of pions (and K and r\ mesons) yields I(m 2 ), 
which, in detail, depends on how xPT is renor- 
malized. (The constants cancel the scheme de- 
pendence.) In a mass-independent scheme 

I(m 2 ) = ?n 2 ln(m 2 //i 2 )- (38) 

All schemes take this form when m 2 <C p?; they 
differ when m 2 ~ p 2 . The natural choice for p is 
500-700 McV, but in most lattice calculations the 
"pion" mass is this large too. The extrapolation 
then depends on how one treats /(m 2 ), and that 
is why chiral extrapolations are uncertain. 

The problem with the chiral extrapolation is 
shown in Fig. 5, which plots £/ vs. r = m™'/m s = 
TOpg/(2mJf — ra 2 .). The black open points come 
from JLQCD's simulations with n/ = 2 [78], re- 
plotted from the original such that = 1 when 
r = 1. These data are all at r > 0.7, and we 
would like reach r = 0.04. The curves are two fit 
Ansatze: a linear fit [78], yielding = 1.13; and 
a chiral- log fit, in which the slope is used to de- 
termine the constant [16], yielding £f — 1.30. 
Both fits are extreme: the linear fit incorrectly 
omits the pion cloud, the chiral-log fit trusts the 
form (38) even when mps > 500 MeV, where 
%PT is a poor description of hadrons [79]. Be- 
cause these data alone do not verify the chiral 
form (38), the uncertainty in £f is large. 
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The gray points in Fig. 5 arc preliminary re- 
sults from HPQCD [80] on the MILC ensembles, 
with ri/ = 2 + 1 and r < 0.8. Although the statis- 
tical errors are still too large to get too excited, 
it is striking how closely they track the chiral-log 
curve (which is fit to JLQCD only). Taking just 
the smallest mass point suggests £/ > 1.2. 

For the benefit of those interested in the CKM 
fit, the rest of this section recommends values for 
the mixing parameters. Let us start with the bag 
factors, whose chiral logs are are multiplied by 
1 — 'ig 2 w —0.05 and, thus, are small: 

Bb. = 1.31 ± 0.10, B Bd = 1.26 ±0.10, (39) 
£ B = 1-022 ± 0.018, (40) 

symmetrizing JLQCD's ranges [78]. These val- 
ues are probably robust, because the bag factors 
seem to be insensitive to nj [57], and the chiral 
uncertainty is a fraction of the total. 

Recommendations for f Bg , /s d , and £ are less 
straightforward. Jb, should be gold-plated [4,5], 
but the comparison of JLQCD and HPQCD 

f Ba = 215(9)(14) MeV (n f = 2) [78], (41) 
f Bs = 260(7) (29) MeV (n/ = 2 + 1) [80], (42) 

is unsettling. The first error bar is statistical; the 
second comes from systematics, such as matching, 
that are partly common. A further uncertainty 
in Eq. (41) comes from converting to MeV with 
the mass of the unstable p meson. Converting 




8Q r i 
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Figure 5. £/ vs. r unquenched. Open symbols 
denote nf = 2 Wilson sea quarks [78]; solid gray 
symbols denote rif = 2 + 1 improved staggered 
sea quarks [80]. The solid line (dashed curve) 
is JLQCD's linear fit (a chiral-log fit [16]) to 
JLQCD's data. Statistical errors only. 



with T splittings would increase f Bs , perhaps by 
15% [81]. Alas, Eq. (42) is still preliminary. It is 
impossible to weigh these issues objectively. My 
recommendation (for this year) is 

f B , = 240 ± 35 MeV. (43) 

To estimate £, the ratio 3 = £/£b fir I fx is use- 
ful. The chiral log largely cancels between £f 
and fir/ fa [82], and the sensitivity to g 2 partly 
cancels between and £ B [16]. From JLQCD's 
data [78] I obtain £ = (/k//tt)H = 1.23 ± 0.06; 
HPQCD quotes (preliminarily) 1.22-1.34 [80], 
though most of the uncertainty here is statisti- 
cal. I shall quote round numbers 

£ = 1.25 ±0.10 (44) 

with the fervent hope that next year's reviewer 
can quote a smaller, yet robust, error bar. 

7. SUMMARY AND PERSPECTIVE 

Most of this review has covered methodology. 
It may be dull, but I hope that this focus will 
help us to meet the challenge of flavor physics. 
The principal concerns — the quenched approxi- 
mation, heavy-quark discretization effects, and 
the chiral extrapolation — are easy enough to list. 
Progress in unquenched calculations [4,5] says 
that we must confront the other two problems. 

Here, and elsewhere [24] , I have argued that we 
should attack both problems by separating scales 
with effective field theories. This gives us a frame- 
work that is theoretically sound, and familiar to 
other theorists, and experimenters too. As long 
as we are not naive in applying effective field the- 
ories, we have good reason to believe that our 
error bars will be robust and persuasive. 

Finally, although it is clear that lattice QCD 
has solid foundations, let us remember that we 
carry out numerical simulations. These are not 
easy for outsiders (even other experts) to grasp 
fully. It therefore always helps to have tests of the 
whole apparatus. In the case of heavy quarks, a 
good set of checks come from quarkonium, where 
the spectrum and also many electromagnetic de- 
cay amplitudes are well measured and (for us) 
gold-plated. Even better for B physics are up- 
coming measurements in D physics. CLEO-c will 
soon measure leptonic and semi-leptonic D de- 
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cays to a few percent. Assuming CKM unitarity, 
one then has a measurement of decay constants 
and form factors. If we can come to grips with the 
chiral extrapolation, lattice QCD has a chance to 
predict the results. Let's not squander it. 

I thank C. Davies, A. Gray, S. Hashimoto, P. 
Lepage, P. Mackenzie, C. Maynard, T. Onogi, S. 
Sharpe, J. Shigemitsu, M. Wingatc, and N. Ya- 
mada for useful discussions. Fcrmilab in operated 
by Universities Research Association Inc., under 
contract with the U.S. Department of Energy. 
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